IPPP/01/53 
DCPT/01/106 
20 February 2002 



Q Prospects for New Physics observations in difFractive 
o ' 

2 : processes at the LHC and Tevatron 

^ '. V.A. Khoze'^, A.D. Martin", and M.G. Ryskin'^'^ 



o 



X 



Department of Physics and Institute for Particle Physics Phenomenology, University of 

Durham, Durham, DHl 3LE 
* Petersburg Nuclear Physics Institute, Gatchina, St. Petersburg, 188300, Russia 



Abstract 

We study the double-difFractive production of various heavy systems (e.g. Higgs, dijet, 
tt and SUSY particles) at LHC and Tevatron collider energies. In each case we compute 
the probability that the rapidity gaps, which occur on either side of the produced system, 
survive the effects of soft rescattering and QCD bremsstrahlung effects. We calculate both 
the luminosity for different production mechanisms, and a wide variety of subprocess cross 
sections. The results allow numerical predictions to be readily made for the cross sections 
of all these processes at the LHC and the Tevatron collider. For example, we predict that 
the cross section for the exclusive double-diffractive production of a 120 GeV Higgs boson 
at the LHC is about 3 fb, and that the QCD background in the bb decay mode is about 4 
times smaller than the Higgs signal if the experimental missing-mass resolution is 1 GeV. 
For completeness we also discuss production via 77 or WW fusion. 



1 Introduction 



Double-diffractive processes of the type 



pp 



X + M + Y, 



(1) 



can significantly extend tlie pliysics programme at fiigli energy proton colliders. Here M rep- 
resents a system of invariant mass M , and the + signs denote the presence of rapidity gaps 
which separate the system M from the products X and Y of proton diffractive dissociation. 
Such processes allow both novel studies of QCD at very high energies and searches for New 
Physics. From a theoretical point of view, hadronic processes containing rapidity gaps play 
a crucial role in determining the asymptotic behaviour of the cross section at high energies. 
From an experimental viewpoint, the presence of rapidity gaps provides a clean environment 
for the identification of a signal. In such events we produce a colour-singlet state M which is 
practically free from soft secondary particles. 

Double-diffractive exclusive processes of the type 



where the protons remain intact, are even better. They allow the reconstruction of the "missing" 
mass M with good resolution, and so provide an ideal way to search for new resonances and 
threshold behaviour phenomena. Moreover, in exclusive processes with forward protons, as 
shown in Fig. 1(a), the incoming gg state satisfies special selection rules, namely it has Jz = 0, 
and positive C and P parity. Hence only a subset of resonant states M can be produced, 
in particular 0"''+ (but not, for example, l"^"*"). Furthermore, the selection rules control the 
threshold behaviour of the production of a pair of heavy particles. For example, the Born cross 
section for the production of a fermion-antifermion pair is proportional to f3^, where f3 is the 
fermion velocity, while for scalar particles the threshold behaviour is just f3. Thus, for instance, 
in this way we can distinguish between scalar quark, q, and gluino, g, pair production. Also, the 
selection rules are of crucial importance in suppressing the bb QCD background when searching 
for the H ^bb signal I, |]. 

We may write the cross sections for processes (|^) and (H) in the factorized form 



where a is the cross section for the hard subprocess which produces the system of mass M, 
and C is the effective luminosity for production at rapidity y. The luminosity C may refer to 
any colourless t-channel states which transfer energy across the rapidity gaps without radiating 
secondary particles. For example, the colourless exchange may be a two-gluon 'hard Pomeron' 
state, or a phenomenological 'soft' Pomeron, or even a. j or W boson. 

The aim of the present paper is to provide a comprehensive way to estimate the numerical 
size of the cross sections for all these double-diffractive mechanisms for producing a wide range 



p + M + p, 



(2) 




(3) 



1 



of heavy systems M (for example, Higgs, dijet, tt and SUSY particles). First, in Section 2 
we determine the luminosity for the different double-diffractive production mechanisms. We 
present plots which readily give the luminosity as a function of the mass M and rapidity y of the 
heavy system, for pp{pp) collider energies ^/s = 2,8 and 14 TeV. Then in Section 3 we give the 
formulae for, and numerical estimates of, the hard subprocess cross section a for the different 
heavy systems M. Armed with this numerical information of C and a, we can immediately 
estimate the observable cross sections cr of a wide variety of different double-diffractive processes 
at the LHC and the Tevatron, see (|]). 

We consider the following double-diffractive mechanisms. Production of the heavy system 

(a) via the exclusive process pp p + M + p, shown in (|]), 

(b) via the inclusive process pp X + M + Y, shown in (|l|), 

(c) via the inelastic collision of two Pomerons. 

The three different mechanisms are shown in Fig. 1. Since we are considering the production 
of a massive system M at high collider energy, the luminosity of the first two mechanisms 
may be calculated from perturbative QCD, using the known parton distributions of the proton. 
Although the cross section a for the hard subprocess is factored off, recall that there may exist 
specific selection rules for the different configurations. In particular in the case of exclusive 
production only the projection with = and P = +1 contributes to a. 



(a) exclusive (b) inclusive (c) inelastic PP 




Figure 1: Different mechanisms for the double-diffractive production of a system of mass M in 
high energy proton-proton collisions. 

We need to be careful in interpreting the cross section a for the subprocess producing the 
massive system, gg — > M. It corresponds to colour singlet production. That is the subprocess 
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amplitude is averaged over the colour indices of the incoming gluons 

= 5Z Sab, (4) 

C a,b 

where Nc = 3 and Aiab is the amplitude for the fusion subprocess g°'g'^ M. All other colour 
factors are included in the luminosity. Moreover, in the case of exclusive double-diffractive 
production, the amplitudes (rather than the cross sections) are averaged over the two transverse 
polarisations of the incoming gluons 

M = It. (5) 

As a consequence, at leading order, we have only Jz = production. On the contrary, for 
inclusive production, the subprocess cross section is averaged over all polarisation states of the 
incoming gluons, as usual. Finally note that for production via soft Pomeron-Pomeron fusion, 
the cross section a includes the convolution with the parton distributions of the Pomeron. 



2 Luminosities for double-diffractive processes 

In this section we collect together the formulae necessary to compute the luminosity functions 
for the processes shown in Fig. 1. We write the effective^ luminosities, for producing a system 
of mass M and rapidity y, in the form 

with i = a,b,c corresponding to the processes shown in Fig. 1. Here we have integrated 
over the transverse momenta of the outgoing protons or outgoing proton dissociated systems, 
according to whether we are considering exclusive or inclusive double diffractive production. 
An important ingredient in the calculation of the luminosity is the inclusion of the survival 
probability of the rapidity gaps to, first, soft rescattering of the interacting protons and, second, 
to QCD radiation. The latter, which results in a Sudakov-like suppression, is included in the 
expressions for L^*^ below. The soft rescattering effects are symbolically denoted by a factor 
S'^ in (^). In practice, we calculate the effects using a two-channel eikonal. As a consequence 
d^) does not have a factorized form, and should be viewed as the soft survival probability 
appropriately averaged over the channels P, The value of the survival factor S"^ depends on 
the particular double-diffractive production mechanism, and may be a function of y and M^. 



^The Pomeron should be regarded as including all multi-Pomeron effects. 
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2.1 Exclusive double-difFractive production 



For the exclusive process shown in Fig. 1(a) we have, to single log accuracy, 



- excl 



TT 



XN'c 



l)b 



(7) 



where h is the t-slope corresponding to the momentum transfer distributions of the colliding 
protons 



dtidt2 



oc e 



where we takeQ 6 = 4 GeV~^. The quantities fg{x,x',Q^, fx'^) are the generalised (skewed) 
unintegrated gluon densities of the protons. The skewed effect arises because the screening 
gluon [Qt) carries a much smaller momentum fraction x' <C x. For small \x — x'\ the skewed 
unintegrated density can be calculated from the conventional integrated gluon g{x,Qf) 0. 
However the full prescription is rather complicated. For this reason it is often convenient to 
use the simplified form p 



fg{x, X , , fx ) 



R„ 



d 



dlnQi 



T{Qt,fi) xg{x,Ql 



(9) 



which holds to 10-20% accuracy. The factor Rg accounts for the single \ogQ^ skewed effect P]. 
It is found to be about 1.2 at the LHC energy, and 1.4 at the Tevatron energy. The Sudakov 
factor T((5t,/i) is the survival probability that a gluon with transverse momentum Qt remains 
untouched in the evolution up to the hard scale = M/2 



T{Qt,fi) = exp 



Qt 



2ti 



l-A 



dz 



(10) 



with A = kt/{fi + kt). The square root arises in (P) because the survival probability is only 
relevant to the hard gluon. It is the presence of this Sudakov factor which makes the integration 
in (^ infrared stable, and perturbative QCD applicable. 

The luminosities M'^dC/dydM^ calculated from (^ and (^ for the exclusive double-diffractive 
production of a system of invariant mass M and rapidity y are shown by the solid continuous 
curves, denoted excl, in Figs. 2, 3 and 4. Fig. 2 shows the dependence of the luminosity on 
M at y = 0, and Figs. 3 and 4 show the dependence on y for produced masses of M = 120 
and M = 500 GeV respectively, at various pp (or pp) collider energies. We use the MRST99 

^If we were to adopt a Regge interpretation, then the 'Pomeron' would be represented in (^ by an uninte- 
grated gluon distribution fg which, for a large hard scale fj,, is described by DGLAP evolution. As a consequence 
of the strong kt ordering, the position of the gluons in impact parameter space is frozen, and hence there is 
no shrinkage of the diffraction cone. This would mean that the corresponding Pomeron trajectory would have 
zero slope, a' = 0. Therefore we choose a constant t-slope, b, which characterises the t dependence of the 
Pomeron-proton vertex. The value 6 = 4 GeV~^ is taken from the fit to the soft hadronic data of Ref. ||^. This 
value is consistent with that observed in J/iJ; diffractive production at HERA 
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Figure 2: The luminosity M'^dC/dydM'^ versus M, for the double-diffractive production of a 

heavy system of mass M with rapidity y = 0. The three plots are for pp [or pp) collider energies 
of = 2, 8 and 14 TeV. Various production mechanisms are studied: the curves marked excL, 
incL, C-inel and soft IP IP correspond, respectively, to production by the exclusive process 
pp ^ p + M + p of Section 2.1, to production by the inclusive process pp ^ X + M + Y of 
Section 2.2, and to production by the processes shown in Figs. 6(a) and 6(b) as described in 
Section 2.3. The 77 luminosity is obtained as described in Section 2.4. 
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Figure 3: The luminosity M'^dC/dydM'^ versus y for the double- diffract ive production by 
various mechanisms of a heavy system of mass M = 120 GeV at ^/s = 2,8 and 14 TeV. The 
notation for the curves is as in Figs. 2 and 6. The upper curve in each plot shows the inelastic 
luminosity (A?7 = 0) assuming that the fusing gg pair are in a colour singlet state. 



partons p to calculate the unintegrated gluon distributions (H), and we calculate the 'soft' 
survival factor 5*^ using the formalism of Ref. P]. We find = 0.045, 0.026 and 0.020 for 
collider energies y/s = 2,8 and 14 TeV respectively. 
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Figure 4: As for Figure 3, but for the production of a system of mass M = 500 GeV at pp 
collider energies of y/s = 8 and 14 TeV. 

The luminosity is presented as the number of effective gluon-gluon collisions per pp inter- 
action. In Figs. 2 and 3 we denote it as gg^^ to indicate that the hard gluons, which interact 
to form the system M, originate within overall colourless (hard Pomeron) t-channel exchanges, 
see Fig. 1(a). This is precisely the quantity which must be multiplied by the cross section a of 
the Jz = 0, colour-singlet hard subprocess gg M, to form the double-diffractive cross section 
a of (I). 

As can be seen from Fig. 2 the luminosity decreases with increasing M, but grows with the 
collider energy ^/s. The reason is that there is an increasing number of gluons as x becomes 
smaller, and because for larger M the double-logarithmic Sudakov suppression (|T0|) is stronger 
due to the higher scale /i = M/2. 

2.2 Inclusive double-diffractive production 

Exclusive production has by far the cleanest signal. Unfortunately, the luminosity, and hence 
the predicted event rate, are small. For this reason we consider inclusive double diffractive 
production, process (|l|), with less restrictive kinematics due to the dissociation of the incoming 
protons, see Fig. 1(b). The cross section can be expressed as exclusive production at the 
parton-parton level 

—* ai + M + a2, (11) 
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Figure 5: The amplitude of Fig. 1(b) multiplied by its complex conjugate, which gives the cross 
section for the inclusive double-diffractive production of a system M . The effective parton 
densities, G{xi) of (|1^), are to be evaluated at scales 



convoluted with the probabilities to find partons ai,a2 in the incoming protons. The process 
is shown in Fig. 5, where the probabilities are denoted by the effective parton densities G{xi). 
At the parton-parton level, ([TT]), the unintegrated gluon distributions in partons oi, 02 may be 
calculated perturbatively in terms of the non-forward BFKL amplitudes Ai,A2. This means 
that we must replace the densities fg in by Ai^/Ti, where = T{kit,fi) is given by (|10D . 
The non-forward BFKL amplitudes Ai are of the form [ill] 

A, = exp(-n,/2)$(F,), (12) 

where ni is the mean number of gluons emitted, with transverse momenta in the range {Qti kit) 
and rapidities in the interval Ar^j 

„..!^A.ln(|). (13) 

The remaining factor ^{Yi) accounts for the usual BFKL (single) logarithms. For rapidity gaps 
Arji < 4 we have 

^ ^ Ar/, < 0.4, (14) 
Ztx 

and, for Ql ^ kf^, it is sufficient to retain only the 0{Yi) term, which gives [|n 



^Yi) ~ 1 + ^ ~ 1.1 ± 0.1. (15) 

We use this value in our numerical predictions. Note that the BFKL amplitudes, which describe 
elastic parton-parton scattering, already contain the double-logarithmic Sudakov suppression 
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which reflects the absence of secondary gluon emission with transverse momenta pt up to the 
momentum transfer ku. The remainder of the hard suppression factor, which accounts for the 
absence of emission in the interval (ku,^), has been incorporated separately in terms of the 
survival factors Tj. 

The dominant (leading log) contribution to the inclusive amplitude of Fig. 1(b) comes from 
the asymmetric configuration -C kf^. Thus the transverse momenta of the two hard active 
gluons can no longer be approximated by ka ~ Qt, as they were in the exclusive case in (0) 
1^. As a consequence each hard gluon propagator, together with its polarization factor, can 
no longer be written as Qt/Qli but must remain ka/kf^. Therefore, the factor l/Q^ in the 
amplitude in (^ now becomes l/{Qtkitk2t)- Moreover, in the hmit <^ kf^ we have ku ~ /c^^. 
Thus, instead of the factor 1/6, which came from the ti integral limited by the proton form 
factor as in (^ and (||), we now obtain the logarithmic dkf^/kf^ integrals. Here we have used 

U = {Q- h)' ^ ~kl ^ -k'l (16) 

Thus, to compute Fig. 5 we must evaluate 

dQl dQf dkl dkj 

Qt Q't ^It 



^ = / ^ 9^ (a,a,a;a^) ^T{r[T2V2. (17) 



In this way, we find the luminosity, (|]), for the inclusive process of Fig. 1(b) and Fig. 5, in 
which the initial protons dissociate and the system M is produced with rapidity gaps A?7i and 
Ar;2 on either side, is, to leading log accuracy, expressed in the form 



i-. = /' G(.,) ^ /' G(..) ^ 4 { I, (18) 

with Nc = 3. The primed quantities arise because the luminosity is obtained by multiplying 
the inclusive amplitude by its complex conjugate, as in Fig. 5. The effective parton densities, 

16 

G{xi) = Xig{xi,kl) + — 51 {(li^i^ ^it) + Q{xi,kl)y (19) 



are integrated from 



xf- = Key + !^ e(^+^''») (20) 



up to 1. These limits ensure that no recoil jets in the dissociated states X and Y lie within the 
rapidity gap intervals Arji and A?72 respectively. 

Due to the asymmetric configurations of the t-channel gluons, Qt <C k^, we have, besides 
Arji, a second logarithm In^kft/Qf) in the BFKL evolution. Using this double log result (that 
is setting $(li) = 1 in (p!2|)), we may simplify expression (0) for X. Then we perform the 
and Q'^ integrations and find 

1 f dti dt2 f 3as . 



t2 



t{\JWIi^t{\Jm,i^, (21) 



where A?7 equals Arji if > |t2|, but equals Ari2 when < |t2|- 

The luminosity M'^dC/dydM^ calculated from ([18|) and @ for inclusive double-diffractive 
production of a system of mass M and rapidity y is shown in Figs. 2 and 3, by the thin 
continuous curves corresponding to two choices of rapidity gaps, namely Arji = Ari2 = 2 and, 
Arji = Ar]2 = 3 in Fig. 1(b). Again we use MRST99 partons 0, to calculate (p^Of). Since we 
are working at the partonic level, a Monte Carlo simulation would be required to obtain the 
precise experimental prediction. Unfortunately, the lower curve for Arj = 3 is more relevant, as 
after hadronization the smaller rapidity gaps with At] = 2 may be hard to identify. 

The effective luminosity, that is the number of gg collisions, is larger for inclusive than ex- 
clusive production due to the larger available phase space, since now the transverse momentum 
of the heavy system M is no longer limited by the proton form factor. Moreover partons of 
larger ku tend to have large momentum fractions Xj, see (|20|). For partons of relatively large Xj 
we sample mainly the smaller size component of the proton wave function, which in turn has 
the smaller absorption cross section 0]. This leads to a larger 'soft' survival probability S"^. To 
be precise we use model II of Ref. |Q to calculate S*^. For inclusive production, the survival 
factor depends on the values of and Ar], since these alter the values of Xi sampled. 



2.3 Production via Pomeron-Pomeron fusion 

The inclusive process studied in the previous subsection allowed secondary particles in the 
proton fragmentation regions, see Fig. 1(b). An alternative possibility is to detect both elasti- 
cally scattered protons, but to permit secondaries in the central rapidity interval 6ri containing 
the heavy system, as schematically sketched in Fig. 1(c). In lowest order perturbative QCD, 
this process is described by Fig. 6(a), which is similar to Fig. 1(a) for exclusive production. 
The only difference is that now we have bremsstrahlung associated with the hard subprocess 
gg^^ M. Therefore we can compute the cross section from the same formula (|^), except 
that the Sudakov-like suppression is weaker. Now the 1 — A upper limit in the z integration in 



([l0|) is specified by 



A = f — cosh ( ^ ) . (22) 




Clearly the luminosity is larger and increases with Srj, as is seen from the dashed curves in 
Figs. 2, 3 and 4. We label the curves^ C-inel. to indicate central-inelastic double-diffractive 
production of the heavy system, see Fig. 6(a). As can be seen from the figures, we do not gain 
much luminosity by allowing extra emission in the central region with Sr] = 2. If we enlarge drj 
up to 4 we enhance the luminosity by an order of magnitude, but at a great price, as will now 
be discussed. 



•^Actually we have included the luminosity of the exclusive process of Fig. 1(a) in the C-inel. results. Strictly 
speaking, therefore, a more precise notation might have been C-total, to indicate that the 'elastic' contribution 
is included. 
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(a) Central-inelastic 



(b) soft PP 



(c) Low-Nussinov 




Figure 6: Double-diffractive production of a heavy system, shown by the bold central arrow, 
accompanied by secondaries in the central region, with forward going protons. 



First, it is important to note that for processes of the type of Fig. 1(c), the mass M 
determined by detection of the forward protons is no longer equal to the mass of the heavy 
system. In this subsection we denote the former mass by Mpp and the latter by Mh- We see 
that Mpp > Mh due to the presence of secondaries in the central region. This clearly applies 
to the process of Fig. 6(a) which we have just discussed, but also to the soft IPIP fusion process 
of Fig. 6(b). For these inelastic IPIP processes we loose all the advantages of the exclusive 
process. The clean central environment of the heavy system is now populated by secondaries. 
The mass can therefore no longer be measured by the missing-mass method. Moreover, 
there is no = selection rule to suppress the QCD background. 

Since the soft Pomeron-Pomeron fusion process of Fig. 6(b) is at present frequently discussed 



12| , |T^], we study it in more detail to confirm that it is indeed not a viable search mechanism. 



In terms of Feynman graphs. Fig. 6(b) corresponds to Fig. 6(c) in which the soft Pomerons have 



been replaced by two-gluon exchange, as originally considered by Low WM and Nussinov |15 



Interestingly, we see that the cross section for the process in Fig. 6(c) is suppressed by a factor a| 
in comparison to that of Fig. 6(a). In the cross section for Fig. 6(b) where, following Ref. |16], 
the Pomerons are treated as real particles with their own parton distributions, the factor is 
reflected in the small values of the effective Pomeron structure functions which are determined 
phenomenologically from diffractive data. It therefore may not be quite so surprising when 
we show below that the soft Pomeron-Pomeron fusion process of Fig. 6(b) has a smaller cross 
section than exclusive double-diffractive production of Fig. 1(a). The cross section^ for soft 



Pomeron-Pomeron fusion was also recently discussed in Ref. ||T2 



*The cross section was also calculated in llSl, where a larger value was found, which wc cannot reproduce. 
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The luminosity factor for the soft Pomeron-Pomeron fusion process, shown in Fig. 6(b), is 

/ O-n \ 2 /• / 1 \ 2("(il)-l) / 1 \ 2(a(t2)-l) 

where we assume the Donnachie-Landshoff parametrization of the elastic pp amplitude |[T^ 

Im A{s,t = 0) = satotis) (24) 
a{t) = 1.08 + 0.25 t 

with s and t in units of GeV^. is the proton form factor and ctq = 21.7 mb. The mass and 
rapidity of the produced system are given by 

r2 1 , f^i 



M'pp = xi X2 s, ^ = 2 [1^2) ' ^^^^ 

For this production mechanism the cross section a in (^, which multiplies the luminosity, 
is the standard inelastic cross section given by the convolution of the parton distributions (of 
the Pomerons) af with the cross section for the aia2 M hard subprocess 



a = 

0-1,0.2 



E / — / — (^1) ^2a!'iz2) cj{a,a2 ~^ M), (26) 



with a = g,q. To estimate a we assume gluon dominance and that zg^{z) ~ 0.7 for our 
kinematic domain, which is consistent with the distribution obtained by the HI collaboration 
from the analysis of their diffractive data [0. Then ( P^ gives 

~ {0.7fa{gg^M){6yy/2 



Zi J Z2 

(27) 



where 

5y = \a{Mlp/Ml). (28) 

Here, to be definite, we have assumed that the heavy system M is a Higgs boson of mass 
Mh- Next, it is important to note that (j{gg H) for the inelastic subprocess within the soft 
IPIP fusion mechanism of Fig. 6(b) is a factor 1/2{N^ — 1) smaller than the corresponding 
(^{99 ~^ H) for the exclusive process of Fig. 1(a). Recall that this factor arises because the 
exclusive process proceeds via the coherent fusion of gluons of different spin and colour, which 
leads to an enhancement of 2(A^^ — 1) in comparison with the inclusive process of Fig. 6(b), 
see (I) and (1). 
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In addition, at Tevatron energies, after observing the leading outgoing proton and antipro- 
ton, the available phase space 5y becomes rather small, 5y ^ 1. In summary, the cross section 
a multiplying the luminosity in is a factor of about 

^ ' 2{Nl-l) 2 - 64 ^^^^ 

smaller for production via soft IP IP fusion (Fig. 6(b)) than via the 'clean' exclusive process 
of Fig. 1(a). Now from Fig. 3, we see that, for the production of a Higgs boson of mass 
Mh = 120 GeV, the soft IPIP luminosity calculated from ( ^31) is only a factor of about 10 
larger than the luminosity for exclusive production. In both processes we use the same survival 
factor]^ 5*^, which was given at the end of subsection 2.1. Thus we obtain our advertised 
result: the expected rate of double-diffractive Higgs production at the Tevatron in the 'dirty' 
environment of inelastic IPIP fusion (Fig. 6(b)) is in fact sma//erQ than the rate for production 
by the 'clean' exclusive process of Fig. 1(a). 

Of course, at LHC energies the available phase space 5y is larger. Recall, from ([27|), a grows 
as {SyY- On the other hand, in going from Tevatron to LHC energies, the ratio of the exclusive 
to soft IPIP luminosities increases by more than an order of magnitude (see Fig. 2), due to the 
growth of gluons at small x. We conclude that there is nothing to be gained by studying the 
soft IPIP fusion mechanism, at least from the point of view of Higgs searches. 

2.4 Production via 77 and WW fusion 

For completeness we give the luminosity factors assuming that the rapidity gaps are due to the 
production of the (M, y) system via 77 and WW fusion. We have 

= ('iV I dh^^^^^^^ {F''^{ti)f f dh^^^^^l^ {F''^{t2))\ (30) 

\Vr/ Jti^nun tl "'i2, min ^2 

where |ti,min| = ^^^m^, m^r is the proton mass and a = 1/137. As the dominant (logarithmically 
enhanced) contribution comes from very small \ti\, that is large impact parameters, here we 
would expect the survival factor S"^ to be close to 1 0. Indeed, we find this to be the case, 
as illustrated by the values of S'^ that are given in Section 3.1.3 for Higgs production via 77 
fusion. The 77 luminosity, (H), shown in Fig. 2 is calculated^] from L^'^, and includes the gap 
survival factor S^. 

^Another way to estimate S*^ is to compare the data with the theoretical prediction which does not account 
for the soft rescattering effects. Special care should be taken if this method is used. If the 'bare' theoretical 
cross section is calculated using, for example, only the mechanism of Fig. 6(b) and does not include the more 
important central inelastic production mechanism of Fig. 6(a), then S'^ = data/theory will be overestimated. 

^This estimate of the soft Pomcron-Pomeron contribution is in reasonable agreement with that of Ref. p^ . 

"^Besides the proton form factor we also include the small contribution from when we evaluate 

F<='"(i) in @. 



13 



The effective luminosity for producing tlie (M, y) system by WW fusion is given by 



■WW 



TT( \ ^ . . dxi /•! dX2 
U{xi) / D{x2) h / D{xi) / U{x2) 

Xi Jx™ X2 Jx5°'" Xi Jx™' X2 



(31) 



w 



wliere = SM^Gp/V^ and Mw is tlie mass of the W boson. The lower limits of the 
integrations are given by (pOl). The effective parton densities, U{xi, kf^) and D{xi, kf^) are 

U{x) = xu + xd + xs + xc 

(32) 

D{x) = xu + xd + xs + xc. 

The survival factor S"^ is calculated using the two-channel eikonal approach described in Ref . 0] . 
We use the more realistic model II [Q, as in Section 2.2. 

Predictions for the WW luminosity are shown in Figs. 3, 4 and 7. Since WW fusion can 
only mediate inclusive production, the missing mass technique cannot be used. However, here, 
the enlargement of the rapidity gaps does not reduce the luminosity as much as it did for gg 
fusion. Fig. 7 shows the WW luminosity for Arj = 2 and Arj = 3, in comparison with that for 
conventional inelastic production via WW fusion for which Arj = 0. These results should be 
compared with the results for gg fusion in Figs. 3 and 4. The upper (dotted) curve shows the 
inelastic luminosity {Arj = 0) assuming that the gg pair are in a colour singlet state. If this 
luminosity is compared with the inclusive luminosity for the production of a M = 120 GeV 
system with rapidity gaps of Arj = 3 on either side, then we loose a factor of 10~^ at LHC 
energies, while the corresponding ratio for production via WW fusion is only 0.02. Production 
by WW fusion is less suppressed with increasing Arj for two reasons. First, since the W bosons 
originate mainly from valence quarks, that is the component of the proton wave function with 
the smaller absorption cross section, the soft survival probability of the rapidity gaps, S"^, is 
larger. Secondly, as the W boson is a colourless object, there is no Sudakov suppression due to 
QCD bremsstrahlung. 



3 Cross sections for the hard subprocesses and for double- 
diffractive production 

Selecting double- diffract ive events means essentially that we are working with a gluon-gluon 
collider^ Indeed, the luminosities presented in Section 2 were given in terms of the numbers 

*For completeness, we also consider production mechanisms wliere the rapidity gaps are associated with 77 
and WW fusion processes. 
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Figure 7: The luminosity M'^dC/dydM'^ versus M for the production of a heavy system at 
y = hy WW fusion, with a rapidity gap of size A77 on either side, at pp or pp coUider energies 
of = 2, 8 and 14 TeV. 

of gg coUisions per proton-proton interaction. As emphasized before, there is an important 
distinction between exclusive and inclusive double-diffractive production of a heavy system M. 
In exclusive events the forward protons select a colour singlet, — incoming gg state. In 
inclusive events the Jz — selection rule is absent. Therefore, to predict the cross section 
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(or event rate at the collider) we have to be careful to ensure that we multiply the relevant 
luminosity by the appropriate gg^^ M subprocess cross section. Recall that the PP super- 
script was to indicate that the gluons are mediating a double-diffractive process (loosely called 
Pomeron-Pomeron production). 

In this section we give the subprocess cross sections da/dfl which are relevant to the double- 
diffractive production of various heavy systems M. The differential form is symbolic. For 
example, for the production of a dijet system, with a rapidity gap on either side, we study 
da/dE^, where Et is the transverse energy of the jets. In this way we can readily predict the 
cross sections for double-diffractive production at various collider energies, using 

^"'^ = (m^ (33) 
dydM^dQ \ dydUy dQ ' ^ ^ 

see (^. The expression in brackets is the luminosity calculated in Section 2, which is to be 
multiphed by the appropriate subprocess cross section, which we enumerate below. We study 
in turn, resonance production (e.g. Higgs, Xb), dijet production, 77 production, ti production, 
the production of SUSY particles (e.g. gg, qq), and various soft phenomena. 

3.1 Resonance production 

The subprocess cross section for the double-diffractive production of a O"*"^ resonance R may 
be written in terms of its two-gluon partial width 

The extra factor of | arises in inclusive double-diffractive production due to the absence of 
the Jz = selection rule in this case. It reflects the spin coherence of the exclusive process 
where the incoming gluon polarisations are correlated, whereas for inclusive production, the 
cross section is averaged over the incoming spins in the usual way, that is 



a""""' ~ \M\', <t'"^' ~ \M\^. (36) 

3.1.1 Higgs production 

In the Born approximation, the width F in ( p^ is just the lowest order R —>■ gg decay width 
Fq. For Higgs production we include the NLO correction to algg H), using 




(37) 

l.bToiH^gg). 
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Note that (|37D is not the complete NLO correction to the H ^ gg partial width, since it does 



not include the contribution from the emission of QCD radiation, see, for example, [0, |20 
For a Higgs boson of mass = 120 GeV, we obtain from and (P7|) 

a-^' ^ 5(1 - ^) l.lpb. (38) 



In the intermediate mass range of the Higgs, a depends weakly on Mh, as Fq ~ M^aKMjj) 
largely cancels the M^^ factor in (0). 

Given the subprocess cross section (|38D and the luminosities of Figs. 2, 3 and 4, we can 
readily predict the cross section for double-diffractive Higgs production at the LHC and the 
Tevatron. To be specific let us take Mh = 120 GeV. Then from Fig. 3 we have 

^ ~ 0.6 X 10-3 (LHC) 

(39) 

~ 0.9 X 10"^ (Tevatron). 

Also we see that for Higgs production Ay ^ 5 (LHC) or 2 (Tevatron) will allow a quick estimate 
of the y integration. Finally multiplying by a of ( P^ ) (which takes care of the integration), 
we obtain 

a{pp^p + H + p) ~ 3fb(LHC) 

(40) 

a{pp ^ p + H + p) ~ 0.2 fb (Tevatron). 

These predictions are consistent with our previous estimates 0. The background to this signal 
for the Higgs boson is discussed in Section 3.2.3. 

3.1.2 Xb production 

Another interesting example is Xb{^^^) resonance production, see, for example, [0, H, |22 
Again a is given by (p^, with |2 



r(x6 - gg) = ro(x6 ^ ^7^7) ( i + 9.8 ^ ) = 550 keV, (4i) 



where we have used the lattice result Fq = 354 keV |24]. This gives 



a-' Jf,) 3.8 nb. (42) 

As an example, consider Xb production at the Tevatron. From Fig. 2 we have 

dC 
dy 

which combined with (^2]) gives 

a{pp p + Xb+p) — 120 pb, (44) 
consistent with the predictions of Ref . . 
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~ 1.6 X 10"^ Ay ~ 2, (43) 

y=0 



3.1.3 Higgs production via 77 and WW fusion 



For completeness we estimate Higgs production by 77 and WW fusion. For 77 fusion, (|3^ ) is 
replaced by 

87r2f(/7^77) 

"'^^-^>= Mi '['-Mi)- (^«) 

The factor of 8, in place of 2, arises because for coloured gluons there was an extra factor 
of (A'^^ — 1)^^ due to colour averaging and the absence of the = selection rule. For 
Mh = 120 GeV we have t{H 77) ~ 7.9 keV, and so 



Hll^H) ~ 5 1 - -— O.lpb. (46) 




The 77 luminosity, integrated over y, is 1.0 x 10 ^ and 1.1 x 10 at the Tevatron and LHC 
respectively, which leads to the cross sections 

(t{pP p + H + p)^^ ~ 0.01 fb (Tevatron) 

(47) 

a{pp ^ p + H + ~ 0.1 fb (LHC). 

Recall that the 77 luminosity includes the gap survival factor S*^, which we calculate to be 
0.75 and 0.9 at Tevatron and LHC energies respectively for the production of a Higgs of mass 
Mh = 120 GeV. The cross sections ( ^71) of the present, more detailed, calculation are lower than 
previous results which were derived in the leading log approximation with 5*^ = 1 ||25|, ^, ^ . 

Note that the strong and electromagnetic contributions to exclusive Higgs production (with 
rapidity gaps) have negligible interference, because they occur at quite different values of the 
impact parameter. Also, it is worth mentioning that the 77 — > if fusion mechanism provides a 
natural lower limit for the Higgs exclusive production rate. 

For Higgs production via WW fusion 



5 1 - 16 nb 

V Ml 



(48) 



for Mh = 120 GeV. The coupling = SM^Gp/V^- 

Of course, WW fusion only mediates inclusive Higgs production. Note that there is almost 
no interference between the WW and gg^^ fusion mechanisms. Recall that the former process 



produces Higgs with high transverse momentum, and as such may give a viable signal [^, 
p^] . Moreover, if we select events with rapidity gaps we suppress the QCD bb background, while 
paying less price for the gaps in WW fusion as compared to those in gg fusion (as discussed at 
the end of Section 2.4). 

18 



3.2 Dijet production 



The Bom cross sections, for colour- singlet production of a dijet system of mass M, are (see, for 



example, E52| 



dt 

dt 

da''"'^ 
dt 

da^""^ 
dt 



{99'' 
{99'' 
{99'' 



99) 



9 7ra| 
4 ~E^' 



99) 



6E^ M2 ^ ' 



9 vr«| / 



M2 



£;|,M2 6 




El 



(49) 
(50) 
(51) 

(52) 



where Et is the transverse energy of the jets and (3 = y 1 — 4m2/M2. We choose the scale of as 
to be M/2 to be consistent with the convention we have used for the luminosities in Section 2. 
For ^ E^ we see, that contrary to 0^+ resonance production (|55|), gg production is a factor 
2 larger in the inclusive as compared to the exclusive process. This is the result of averaging 
over the polarisations of the incoming gluons, (|^). Moreover we see that the Born cross section 
(|50|) for exclusive gg^^ — > qq production vanishes as — > [^, which is a consequence 
of the Jz = selection rule for the gluon polarizations 0, |T^. 



Rather than working in terms of the subprocess variables and t, it is often convenient 
to use E^ and Srjj, where Srjj = rji — ri2 is the difference in rapidities of the two outgoing jets. 
It turns out that factors in the Jacobian cancel, so that we obtain the relatively simple relation 
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da 
It I 



'dM^ 



da 
dE^/ 



d{6r]j] 



(53) 



Thus the product of the luminosity (shown in Figs. 2, 3 and 4) and da/dt, will give an evaluation 
of the differential cross section da / dyd{Srjj)dEl. 

The large cross section for exclusive dijet production provides the first opportunity of ex- 
periment to check the double-diffractive predictions, already at the Tevatron. For example, let 
us calculate the event rate for jets in the domain defined by the intervals 25 < Et < 35 GeV 
and Srjj < 2. If we integrate ( |15| ) over these intervals, and multiply by the luminosity dC/dy ~ 
3 X 10"^ appropriate to M = 2i?T cosh((5?7j/2) ~ 70 GeV, then we obtain 



Aa{pp — ^ p + jj + p) ~ 60 pb 



(54) 



corresponding to jets with 25 < Ej- < 35 GeV and \t]i — ri2\ < 2, which is in qualitative 
agreement with previous results fT^, |^. Thus even for a Tevatron luminosity of 1 fb~^, we 
estimate there should be 60,000 exclusive events in this particular domain. 
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3.2.1 Dijets as a luminosity monitor 



There are, of course, uncertainties in our estimates of the cross sections for double-diffractive 
processes. First, the calculation of the soft survival factor S"^ is based on a simplified two- 
channel eikonal model 0, ^ . With the presently available soft diffractive data, these estimates 
of 5^ are as good as possible, but it is hard to guarantee the precision of predictions which rely 
on soft physics. Secondly recall, from Section 2, that the luminosity was calculated in terms 
of the unintegrated, skewed gluon distribution fg, and that C oc fg. Thus a 20% uncertainty 
in fg results in a factor of 2 uncertainty in the luminosity. We calculated fg in terms of the 
conventional gluon distribution obtained from the results of a global parton analysis of deep 
inelastic and related hard scattering data However the prescription is only justified at 
leading order and for small x [§, 0| . In the present paper we use an improved formula for fg , 
and as a result the double-diffractive luminosity at Tevatron energies become a factor of two 
larger than in Ref. 00. However the same luminosity determines the cross sections for different 
double-diffractive processes. Thus dijet production, where the cross section is largest, offers an 
excellent way to monitor the gg^^ luminosity ||^. 

3.2.2 Dijets as a gluon factory 

Another application of the large rate of exclusive dijet production is as a gluon factory. The 
remarkable purity of the colour-singlet, Jz = 0, di-gluon system provides a unique environment 
to study high energy gluon jets 0. 

3.2.3 Dijets as a background to the Higgs signal 

The large dijet production rate may produce a huge background to the Higgs H ^ bb signal. 
Recall that the cross section for exclusive double-diffractive Higgs production was calculated in 
Section 3.1.1. 

In fact, for any 0^+ resonance the ratio of the exclusive R gg signal to the gg background 

is 

where AM is the experimental missing- mass resolution, the scale of as is M/2 and the factor 
9.7 comes from imposing a cut 60° < 6* < 120° in order to improve the signal-to-background 
ratio. 6* is the decay angle in the dijet rest frame. Note that the major part of NLO effects 
and uncertainties in the luminosity cancel in ratio ([551). 

^The unintegrated skewed gluon distribution fg{x, x' , Ql,y?) has been tested, both by providing a successful 
description of the data on diffractive vector meson production at HERA |Q and by agreement with the upper 
bound on diffractive dijet production at the Tevatron However in both cases a much lower scale is 

probed and we did not need precision better than 20-30%. 
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For H — s> 66, the signal-to-QCD background ratio is 

S(gg^^^bb) ^ St. ... f M \^ /250 MeV\ 

^JL_L ^ 4.3 X 10- Br(H 66) (^^) (^^) , (56) 

which is very small 0. Fortunately, if we tag the b and b jets to reject the gg events, we can 
strongly suppress the QCD background. Recall that the gg^^ — > bb QCD background process 
is suppressed by colour and spin factors, and by the = selection rule. In this way the 
background is suppressed by an extra factor 

However the full suppression is only true in the Born approximation. Recall that large angle 



gluon radiation in the final state violates the Jz = selection rule [^, so in the bb QCD 
background the suppression factor ml/Mfj is replaced by 05/71. The final result for the H ^ bb 
signal-to-background ratio is therefore 

Sjgg^^^H^bb) ^ / 250 MeV \ 

B{ggPP ^ bb) - [ AM ) ' ^^^^ 

for a Higgs boson of mass Mh = 120 GeV. 



3.3 Double-difFractive 77 production 

At first sight, the subprocess gg^^ —>■ 77 appears attractive to serve as an alternative gg^^ 
luminosity monitor for the exclusive double-diffractive processes[^. However it turns out that 
the event rate is too small. Using the known results^ for the QED 77 — > 77 helicity amplitudes 



371 , we calculate the Jz = 0, parity P = +1 subprocess cross section to be 

<j(30° <e;< 150°) ~ 0.3(0.04) pb (59) 

for M^^ ~ 50(120) GeV. When multiplied by the luminosity at the Tevatron of 1.8 x 10^^(2 x 
10"'^), we obtain, assuming / dM'^/M'^ ~ 1, 

a(pp^p + 77 + p) ~ 0.5 fb (0.008 fb). (60) 

On the other hand, for the LHC for M^y ^ 120 GeV we expect 

cr(pp — > p + 77 + p) ~ 0.12 fb. (61) 

Note that for M^^ ~ 120 GeV the observable rate of Standard Model exclusive 77 events is 
expected to be much less than even Higgs bb signal. In comparison to exclusive gg^^ gg 
subprocess, the gg^^ —>■ 77 rate is smaller by a factor of 10^ for Etq — Et-^- 

^"We are grateful to Mike Albrow and Beate Heinemann for discussions on this proposal. 
^^We thank Andrei Shuvaev for checking various formulae for 77 production. 
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This obvious disadvantage of the Standard Model exclusive diphoton production can be 
turned into an attractive advantage for searches at hadronic colliders of signs of the existence 
of extra dimensions, which for instance appear in theories of low-scale gravity The small 
background rate, induced by the conventional box diagram contribution, may allow a very 
sensitive high-mass diphoton probe of the effective scale of quantum gravity effects, Ms (see. 



for example, [39, ffl). To illustrate the level of the possible signal we follow the estimates of 



Ref. ||39| with effective scale Mg = 1.5 TeV. At M^^ ~ 500 GeV one then would expect, for 



30° <9*< 150°, 

^^xtL dim. ~ / ) 50 fb, (62) 



where the factor F is given by (see, for example, p9| ) 




F = 2/(n-2) for n > 2, (63) 

where n is the number of extra (compactified) dimensions. We consider here the inclusive 
configuration of Fig. 1(b), since for exclusive LO kinematics there is no = point-like 
gluon-gluon coupling to a "graviton" . 

Taking / dM'^/M'^ ~ 1, after multiplying by the double-diffractive inclusive luminosity of 
2 X 10^^ (for Arj = 3), we would expect at the LHC 

Aaextra dim. {pp X + + Y) ~ 0. IF^ fb, (64) 

which, assuming that F"^ is of order 1, should yield about 10 events, for an integrated luminosity 
oi C = 100 fb~^. It is quite possible that future Tevatron results may increase the existing lower 
limit on Ms, and hence make a diphoton signature for extra dimensions invisible. 

Of course, in the inclusive configuration, there may be some additional background induced 
by 77 emission off the quark lines in the process where the rapidity gaps are created by colourless 
qq t-channel exchange. However for gaps with At] = 3, a preliminary estimatej^ shows that 
this quark contribution does not exceed the background contribution coming from the gg —>■ 77 
box diagram. Thus the total background should not be greater than about 5-10%. 

3.4 ti production 

The exclusive double-diffractive production of tt pairs may provide new opportunities for study- 
ing top quark physics. In particular, it offers a novel probe of the QCD dynamics in the tt 
system. The Born cross section is given by (0). Due to parity conservation, the ti pair is 
produced in a P-wave state, and so the threshold behaviour goes as where /? is the quark 
velocity in the tt rest frame. Of course, the Born cross section is modified by Coulomb and 



top-quark width effects ||41|, ^], but for illustration it is sufficient to neglect these, see our 
estimates at the end of subsection 3.5.2. 

Note that the observation of tt exclusive production can be used as a template of the ability 
of the missing mass method to find signals of New Physics. 
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We plan to study qq exchange in more detail in the future. 
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3.5 SUSY particle production 



Double-diffractive exclusive processes, in principle, provide a unique opportunity to investigate 
the whole of the strong interaction sector of physics beyond the Standard Model. As an example, 
we consider super symmetry, which is a front-runner in searches for New Physics. Here we ignore 
the current theoretical SUSY prejudices and assume for illustration, first, that the gluino is the 
lightest supersymmetric particle (LSP). The second scenario that we consider is that unstable 
gluinos and squarks exist with masses close to the current experimental limits. Note that it 
is difficult to separate gluinos and squarks produced in standard inelastic hadronic collisions. 
Moreover, gluino studies appear especially promising in gg collisions since these cannot be easily 
achieved at linear e~^e~ colliders. 



3.5.1 Gluinoball production 



In recent publications |^ an interesting possibility is discussed that a gluino g is the LSP (or 
next-to-LSP) or, to be specific, that it does not decay within the detector. The allowed mass 
window is claimed to be 25-35 GeV. Within such a scenario there should exist a spectrum 
of gg bound states (gluinoballs or gluinonia), which may reveal themselves in gluon-gluon 
collisions. The orbital angular momentum L must be odd in such gg^^ collisions, due to parity 
conservation. Therefore the lowest lying colourless bound states is the 0"^+('^Po) state, which we 
denote G, with principal quantum number n = 2. Relative to the 2mg threshold, the energies 
of such P-wave gg bound states are 



9 

- m, 
4 



9 9 



(65) 



with n > 2, where as is to be evaluated at the Coulombic scale kc = 3as'mg/2 (see, for 
example, ^). We may estimate the partial width of the G gg decay using the Coulombic 



approximation. The old results of the two-photon decay of P-wave positronium |^ may be 
applied (see, for example, [0, ^), which gives 



T{G^gg) = 6A aUrrig) aUkc) 

0.2 MeV, 



(66) 



.60 GeV 

where Mq ^ 2mg. Then, on using (0), we obtain 



— 18 pb. 



M. 



(67) 



The anticipated cross section for the signal is quite sizeable. For example, for Mq = 60 GeV, 

a{pp^p + G + p) ~ 0.4 pb (LHC), 

(68) 

a{pp —^p + G + p) ^ 20 fb (Tevatron). 
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However the signal-to-background ratio is 

S{ggPP ^G^gg) ^ ^ , /250 MeV\ / 



0-6 X 10-M — -— (69) 



Bigg'''' gg) ^ \ AM J Uo GeV 

which makes detection difficuh, even with angular cuts. Some words of caution are in order 
here. Our estimates are based simply on the lowest-order Coulombic formula, ([661). More 
refined calculations of r(G gg) are certainly needed. Secondly, we have neglected the 
possible interference between the signal and the background. 

The threshold production of (quasi) stable gluino pairs may be strongly affected by QCD 
final-state interactions, in analogy to the celebrated Coulomb threshold phenomena in QED 
i7|. As follows from the results of Refs. ^], in the zero-width approximation the P-wave 



threshold cross section, dap{gg; M), is 

dapCg~9;M) = daP{~g~9;M) \ij-g{0)\' (^1 + , (70) 

where dap\gg, M) denotes the Born cross section at cm energy M > 2mg. "ipgiO) is the Coulom- 



bic wave function of a colour-singlet gg S'-wave state evaluated at the origin (c.f. Refs. HI)- 



with 

Z. = (72) 



where jSg = y 1 — 4m|/M2 is the velocity of the g. It looks natural to choose the scale of as 
in (|72D to be pc = max 

Since Zg > 1, QCD Coulombic effects drastically modify the gg excitation curve, strongly 
enhancing the production rate at or near threshold. According to ([7^), for < An"^ the 
threshold cross section rises with increasing (3g approximately as rather than /5| as one 
would expect from the Born result. Moreover, for Z| 3> 47r^, the colour Coulombic exchanges 
in the final state completely compensate the kinematical /3| threshold factor contained in the 



Born cross section 42 



for ~ 8 GeV which corresponds to AM ^ 1 GeV. This is a dramatic effect relative to the 
lowest-order expectation. The (LSP) gluino-pair production process would lead to distinctive 
signatures in hadronic collisions, in particular to events with jets accompanied by missing pt, 



see, for details, Ref. ^ 



Despite all the attractive features of the gluino-LSP scenario [Q, it is possible that direct 
experimental studies at the Tevatron will close this window. 
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3.5.2 Gluino and squark production 



A more plausible scenario is to search for massive unstable gluinos g or squarks g, which will 
reveal themselves as multijet final states with missing energy and/or leptonsQ. As noted above 
the advantage of exclusive double-diffractive studies is that gluinos and squarks can be separated 
by their respective and /? threshold behaviour. The Born subprocess cross sections are 



da' 



cxcl 



dt 



27 1 7m| m| 
2 2 M2 ^' 



3' 



(74) 
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2m|^ 



(76) 



(77) 



Note that ( |71[ ) and (^) are the same as (|5D|) and ( p^ for colour-singlet gg qq, except for 
the colour-factor 27/2, and a factor 1/2 which refiects the Majorana nature of the two gluinos. 
The (3"^ behaviour of exclusive gluino-pair production arises from (3"^ of the last factor of ([7^) 
and (5 from the integration over t. The cross sections shown in (^) and ([77|) are written for 
one fiavour and one type (L or R) of squark q. 

Near threshold the formula have to be modified to account for sparticle widths and the QCD 
Coulomb interaction 0. Unfortunately the cross sections are very small near threshold. 
For example, if we take sparticle masses of 250 GeV and integrate from threshold (500 GeV) 
up to 625 GeV, then we find 



^^cxci(~~) ^ 6.5 pb. 



Aa""='(g^) ^ 1.8 pb. 



(78) 



For comparison the tt cross section, integrated from threshold up to 437 GeV, corresponding 
to the same value of /3 = 0.6, is 

A<T""'=^(tt) ~ 2.2 pb. (79) 

After multiplying by the double-diffractive luminosities of 2 x 10~^ for gg and 4 x 10~^ for ti 
we have 



Aa{pp ^ p + gg + p) ^ 0.15 fb 



(80) 



Aa{pp ^ p + qq + p) 



0.04 fb 



(81) 



Aa{pp ^p + tt + p) ~ 0.1 fb (82) 

^^At the moment we have no clear understanding which of these sparticles is the hghtest, but the conventional 
belief is that the stop t is the lightest squark. 
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for 'near' threshold production at the LHC It is quite plausible that the masses of the light 
squark flavours are nearly degenerate, as are the masses of and q^. This may allow a higher 
rate of qq events. 

Another possibility to increase the yield of SUSY production is to consider the inclusive 
configuration of Fig. 1(b). Then the effective luminosity (~ 2 x 10"'^) is much larger. The 
corresponding subprocess cross sections, integrated from threshold (500 GeV) up to 625 GeV, 
are 

/\a^''\g~g) ^ "24. pb, Aa'"^='(gt) ^ 1 pb, (83) 

leading to 

Aa{pp X + gg + Y) ~ 50 fb, (84) 
Aa{pp ^ X + q~q + Y) ~ 2 fb. (85) 



3.6 Soft Phenomena 

Finally, we discuss the possibility to study soft strong interactions in diffractive events with two 
(or more) rapidity gaps (see, for example, [5^, |5^). The (total) cross section of the Pomeron- 



Pomeron interaction is proportional to the square of the triple Pomeron coupling g^jp. If we 
take the parameters from the analysis of Ref. then we obtain the subprocess cross section 

atilPIP) ~ 1.2 mb. (86) 

When we multiply by the "soft IP IP" luminosity of Section 2.3, shown in Fig. 2, we predict a 
two rapidity gap cross section 

(72 ~ 1-10 fib, (87) 

depending on the gap size. For the LHC it gives a very large event rate, which may be measured 
in a low luminosity (high j3 optics) run. 

Such events offer an excellent opportunity to answer the many outstanding, interesting 
questions which remain from the intensive studies of soft physics of 25 or more years ago (see. 



for example, |5T|, |5^). A key problem is that the cross section of events containing n rapidity 
gaps is of the form 

a„ ~ ((73iplns)27(2n)!, (88) 

which grows with energy. On summing, we have a cross section a ~ (s)^^-^, which appears to 
violate unitarity. We must find a mechanism to suppress this huge cross section. We recall the 
scenarios which were considered long ago. 

(i) To obtain a self consistent asymptotic theory it was proposed that the triple Pomeron 
coupling vanishes when the momentum transferred through the Pomeron goes to zero, 

gsAt) ^0 as t 0. (89) 
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The diffractive data do not confirm this proposal, although in the presence of strong 
iP-cuts it is still an open question. The study of double-Pomeron-exchange processes in 
low luminosity LHC runs, where the small t region becomes accessible, would clarify the 
situation (see and references therein). 

(ii) It was suggested that the soft survival probability S"^ becomes much smaller for events 
with many gaps. The available CDF data [^, on double and single diffractive production, 
indicate otherwise, but again it should be checked at LHC energies. 

(iii) It is more natural to expect 5^ does not depend strongly on the number of gaps, but rather 
to propose that S"^ decreases rapidly with increasing energy. In this scenario, rapidity gap 
events only occur in peripheral (large impact parameter) collisions where the opacity fl{b) 
is still small. It follows that mean transverse momentum of the secondaries (kt) created 
in the central IP IP collisions of Fig. 1(c) will be smaller than (kf) for the usual inelastic 
interaction at an equivalent energy, ^/Si^l = M [^. This prediction could be checked at 
the LHC. 

We emphasize that processes with 3 or more rapidity gaps have not yet been seen. It is 
clearly crucial to detect such events. To estimate o"3,(T4, . . ., we assume factorization, and the 
same survival factor S*^ for cxi, (T2, era, . . .. Then the probability to observe an additional rapidity 
gap will be 

where we cancel the factor 5*^ in the empirical double-diffractive cross section ace, since the 
survival factor 5*^ is the same for all cr.„. Of course, we have to account for the phase space 
available for each gap. To make a phenomenological estimate, we start with the two-gap process 
of Fig. 1(c) and evaluate the probability to observe a third gap at the Pomeron-Pomeron energy 
^Jspp = M. For the LHC it corresponds to ^Jsjpp ~ 300 GeV. Using the existing data in the 
ISR to Tevatron energy range, together with the S"^ calculations of Ref. [0, we find 

cugap ^ 0.2-0.4 (91) 

and 

c^s - ~ 0.3-3 iih. (92) 

Another possible soft phenomenon which is worth investigation is the "elastic" Pomeron- 
Pomeron scattering process 

JPIP ^ X^Y^, (93) 

where X*^ and Y*^ are glueball states. We call it elastic as the glueballs may lie on the Pomeron 
trajectory. That is there may exist a glueball dominance model of the Pomeron, in analogy to 
the vector meson dominance model of the photon. 
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4 Conclusions 



In this paper we summarize, and extend, the results of our previous studies of double- diffr active 
processes with a rapidity gap on either side of a centrally produced heavy system. We have 
attempted to present the results in a unified, and user- friendly, form so that the expected cross 
sections at the Tevatron and the LHC can be readily estimated, simply by multiplying the 
appropriate luminosity of Section 2 with the relevant subprocess cross section a of Section 3. 
We hope that it provides a useful framework to study double-diffractive phenomena, not only 
in their own right but also as a tool to probe physics beyond the Standard Model. Although 
many cross sections are predicted to be small, there are several interesting processes with viable 
signals. In particular, the exclusive double-diffractive production of an intermediate mass Higgs 
boson at the LHC, and the possibility of using the Tevatron and the LHC as a 'pure' gluon 
factory. For example, for the production of a Higgs boson of mass Mn — 120 GeV at the LHC 
we predict a cross section 



where AM is the experimental missing-mass resolution. 
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^ p -|- i7 -|- p) 
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(94) 



with a signal-to-QCD background ratio of 




(95) 
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